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Many-body entangled systems, in particular topologically ordered spin systems proposed as re- 
sources for quantum information processing tasks, often involve highly non-local interaction terms. 
While one may approximate such systems through two-body interactions perturbatively, these ap- 
proaches have a number of drawbacks in practice. Here, we propose a scheme to simulate many-body 
spin Hamiltonians with two-body Hamiltonians non-perturbatively. Unlike previous approaches, our 
Hamiltonians are not only exactly solvable with exact ground state degeneracy, but also support 
completely localized quasi-particle excitations, which are ideal for quantum information processing 
tasks. Our construction is limited to simulating the toric code and quantum double models, but 
generalizations to other non-local spin Hamiltonians may be possible. 
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Many-body entanglement arising in strongly correlated 
systems is a very promising resource for realizing various 
ideas in quantum information, such as quantum commu- 
nication and quantum computation. In particular, topo- 
logically ordered spin systems can be employed for re- 
liable storage of quantum information inside the degen- 
erate ground space [l| and for fault-tolerant quantum 
computation with non-abelian anyonic excitations |2J. 
These topological approaches may resolve many prob- 
lems in quantum information science; qubits are encoded 
in many-body entangled states and are thus naturally 
protected from decoherence. 

Unfortunately, topologically ordered spin systems ca- 
pable of quantum information processing are very diffi- 
cult to realize physically. Many proposed topologically 
ordered spin Hamiltonians, such as the toric code, quan- 
tum double model [2[, and string- net model [3j], involve 
highly non-local interaction terms; this is a stark contrast 
to Hamiltonians which occur in nature, which have only 
geometrically local two-body interactions. Moreover, the 
resource systems above are known not to be supported 
by any two-body Hamiltonian [J]. 

Many efforts have been made to construct two-body 
Hamiltonians which "approximate" non-local resource 
Hamiltonians. The most commonly used approach is to 
approximate target Hamiltonians through so-called "per- 
turbative gadgets" joHlOj . The central idea of perturba- 
tive gadgets is to design a two-body Hamiltonian whose 
leading pcrturbative contribution gives rise to the de- 
sired many-body Hamiltonian; unfortunately, most ob- 
tained two-body Hamiltonians are not exactly solvable, 
and their properties are hard to determine except for a 
few exactly solvable examples [ll|, [l2|]- In addition, the 
perturbativc Hamiltonian only approximates the target 
Hamiltonian, and may give a very weak effective Hamil- 
tonian with a rather small energy gap. Furthermore, 
quasi-particle excitations (energy eigenstates) arising in 
perturbative Hamiltonians cannot be created through 
completely localized manipulations of spins; excitations 
are always delocalizcd and the ground state degeneracy 



might be split for finite system sizes, resulting in fatal 
errors in practice. While a non-perturbative approach 
based on the PEPS formalism was developed for simu- 
lating the cluster state for measurement-based quantum 
computation [13| , such an approach may not be applica- 
ble to degenerate systems with topological order. 

Here, we propose a scheme to simulate topologically or- 
dered Hamiltonians through two-body interactions non- 
perturbatively. Our scheme builds on previously estab- 
lished ideas in perturbativc gadget studies, such as the 
use of hopping particles proposed by Konig [8| , and the 
encoding of single particles into multiple particles used by 
Brell et al [l(|. Combining these remarkable insights, we 
are able to construct the first topologically ordered spin 
system which satisfies the following; 1) The Hamiltonian 
has at most two-body, geometrically local interactions. 
2) The Hamiltonian has exactly solvable ground states 
and low-energy excitations, and is provably gapped for 
all system sizes. 3) The ground space of the Hamiltonian 
is exactly connected to that of the target Hamiltonian 
through local unitary transformations, and anyonic exci- 
tations are completely localized. 

For clarity of presentation, we illustrate the gadget 
construction for the toric code. A generalization to the 
quantum double model is also possible (see appendix [Cj. 

Modified toric code — We begin by defining a modified 
version of the toric code, also known as the Zi lattice 
gauge model, that we will simulate through a two-body 
Hamiltonian. Consider a system of qubits defined on 
edges of a square lattice with periodic boundary condi- 
tions. Unlike the conventional toric code, two qubits live 
on each edge in our construction (see Fig.[TJa)), governed 
by the following Hamiltonian 
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where s, p and e represent "star" , "plaquette" and "edge" 



respectively, as defined in Fig. [T{b)(c)(d). ^j an d %j 
are Pauli X and Z operators on qubit j, and J is some 
positive constant. The model is exactly solvable since 
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FIG. 1. (a) Construction of the modified toric code. Dots 
represent qubits. (b) A star term A a (red online), (c) A 
plaquette term B p (blue online), (d) An edge term C e (green 
online), (e) Two pairs of logical operators. 

interaction terms commute with each other, and it can 
be considered to be a stabilizer code. The model has 
four degenerate ground states, as in the toric code. In- 
side of the ground space, A s = B p = C e = 1, meaning 
A s \ip) = B p \ip) = C e \ijj) = \ip) for all s, p and e when \ijj) 
is a ground state. Notice that one can create the toric 
code from this model by applying controlled-NOT gates 
between pairs of qubits on each edge. Since the toric code 
and the modified model are connected through local uni- 
tary transformations, the y a re considered to be in the 
same quantum phase [lj, Qjf . The ground space of the 
modified toric code has a four-fold degeneracy, as seen by 
writing down two pairs of "logical operators" which com- 
mute with the Hamiltonian but anti-commute with each 
other (see Fig. [He)). The non- locality of logical opera- 
tors makes the model of great interest for robust storage 
of quantum information. 

As a first step towards obtaining a two-body Hamil- 
tonian simulating this modified toric code, we group the 
four qubits in each plaquette into a single composite par- 
ticle with a 16-dimensional space. While B p becomes 
one-body, and C e is two-body through this grouping, the 
star term A s is only reduced to four-body Below, we 
provide a scheme to simulate A s through only two-body 
terms. 

Gadget Hamiltonian- The central idea behind our con- 
struction is to add a "gadget particle" at each star (see 
Fig-IH a ))- The gadget particle has four possible spin val- 
ues m s = 0, 1,2,3. We replace the four-body star term 



A s with two-body terms Hh op and H S Meid which involve 
the gadget particles: 



Hgadget — H p + H e + Hhop + Hshieid 

H p = — J 2_^ B p , H e = —J 2_^ Ce ■ 

p e 

The hopping term is Hhop = ^2 S h s where 
h s = -U\m s = 0)(m, = 0| - t (Uj + D.) 



(1) 



Dl = S~] \m s + l)(m s \ ®A s (m s ) (mod 4), 

m,=0,l,2,3 



where U and t are some positive constants, and m s rep- 
resents the spin value of the gadget particle at s. Terms 
A s (ra) are products of two Pauli X operators as depicted 
in Fig. Hlb). Since A s {rn) are one-body operators when 
qubits in a plaquette are viewed as a composite particle, 
hopping terms are two-body. This hopping term will ef- 
fectively induce star terms A s since A s = (D\) A . The 
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FIG. 2. Construction of the hopping term Dj. (a) Gadget 
particles at stars, (b) Terms A s {m) that are tensor products 
of two Pauli X operators. Each term acts on two qubits (red 
online), depending on spin values of gadget particles. 

shielding term H s hi e id consists of two-body interactions 
between gadget particles: 

H S Meid = J^T e (ms)T r (m s+ x) + T d (m s )T u (m s+ y) 



where s + x and s + y are unit translations of a "star" s 
in the horizontal and vertical directions, and 

T e (m) = 1 - 2(5 m , 2 , T r (m) = 26 mt o - 1 
T d (m) = 1 - 2(5 m . 3 , T u (m) = 1 - 2<5 mj3 

with Sn-^m' = 1 for m = m! and otherwise. As we will 
see below, this choice of the shielding term decouples ef- 
fective interactions between neighboring gadget particles, 
and makes the model exactly solvable. 

Decomposition into subspaces — Now, we solve the gad- 
get Hamiltonian in Eq. (jDH) . It is convenient to decom- 
pose the entire Hilbert space into subspaces. Let us de- 
note computational basis states whose gadget values are 



all |0)s: \ip(d)) = \0) ga dget ® \d) qu ut where \d) qubit repre- 
sents spin values \dj) for qubits, and \0) qu bit means that 
all the qubits are |0)s. We define the subspace A4(d) such 
that it is spanned by all the states which can be reached 
from \4>{d)) by applying D\: 

M{d) = { Y[(Dl) x ° |V(d)), for all A s >. 



We can verify that M.{d) is an invariant subspace of 
H gadget- Then, one can solve the gadget Hamiltonian 
inside each subspace M{d) independently. 

Ground state subspace — We will first solve for the 
ground state inside -M(0), and then will show its lowest 
energy state to be a ground state. We note that inside 
M.(Q) B p = 1, and thus plaquette terms need not be 
considered. Denoting the total number of stars as N, we 
may view M(0) as the Hilbert space of N particles. 



write the ground state as 
7 



|A) 



iA S >=n(^) As i^(o)>- 



(2) 



Noting that (Dj) 4 = A s , (Dj) 8 = I, these particles can 
be considered to have eight-dimensional Hilbert spaces, 
\ s = 0, . . . , 7 [la]. In this "A- representation" , the hop- 
ping term Hh op can be written as a one-body Hamilto- 
nian: Hhop = y h s where 

h s = - U(\\ s = 0)(A S = 0| + | A. = 4)(A S = 4|) 

7 

-t^2 (\X s + l)(X s \ + h.c) (mod 8). 

A s =0 

However, edge terms C e are not one-body inside M(0). 



A key idea behind our gadget arises from the fact that 
these two-body interactions arising from C' e can be ex- 
actly cancelled by adding the shielding term H s hi e id- In- 
side -M(0), edge terms have the same action as the fol- 
lowing two-body terms involving gadget particles: C e = 
Te(rn s )T r (m s +x) for a horizontal edge e connecting s and 
s + x, and C e = Td(m s )T u (m s +y) for a vertical edge e 
connecting s and s + ij, as one can verify from direct cal- 
culations (see appendix [3}. Then, the edge terms are 
exactly cancelled: H e + H s hi e id — inside AA(Q). This 
means the gadget Hamiltonian is one-body in the "A- 
representation" : H gadget = const + ^ s h a . 

Because of this, all energy eigenstates inside .M(O) can 
be written in the tensor product form \ct) = (^) s \u s ) 
where |a s ) = J2\ a s(A s )|A s ). The lowest energy state is 

IV'Gs(O)) = ® s |a ), where a (A) = a (A + 4) for all A. 
Therefore, returning from the A-representation, we can 



lVW0)> = nE a °( A )(^) A hK0)> 

s A=0 

=n( /+ ^)E a °(Ax^)>(o))- 
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We see that there is a finite energy gap inside M. (0) , since 
H gadget acts as a one-body Hamiltonian. 

Unitary Connection — This lowest energy state 
|'0Gs(O)) is connected to the ground state of the mod- 
ified toric code through the following local unitary 
transformation: 

3 
U = ]JU S , U S =J2 \ m s)(m s \ {J A 8 (m). (3) 

s m s — m<m s 

In particular, we have U\if) GS (d)) = \do)f^dget ® 

\lpTorie(0))qubit where |o7 ) = J2m=0 a o( m ) \ m ) I and 

H'Toric{0))qubit = Ylsi 1 + A s)\0) is a ground state of 
the modified toric code. We may verify that the gad- 
get Hamiltonian has three other ground states \ipas(di)), 
i = 1,2,3, inside M.(di), connected in the same way to 
the ground states |VToric(d»)) of the modified toric code. 

It is then simple to find the logical operators for the 
gadget Hamiltonian; they are those of the modified toric 
code conjugated by U: WXiU, W X 2 U , U^Z-JJ and 
XJ^Z-iU . The ground space is topologically ordered since 
it meets the criteria for the stability against local pertur- 
bations proposed in [l7| ■ 

Anyonic excitations, which are also energy eigenstates, 
can be created by applying "segments" of logical opera- 
tors combined with local operations on gadget particles 
in a similar way to the conventional toric code. As a re- 
sult, excitations can be created only through completely 
localized manipulations of spins in small regions. This is 
in striking contrast to pcrturbative Hamiltonians where 
anyonic excitations arc delocalized, and cannot be created 
through completely localized manipulations of spins. 

Energy gap — Finally, we show that \ipGs(di)} are the 
ground states of the gadget Hamiltonian. To do so, we 
prove that the lowest energy states within other non- 
ground-state subspaces M(d) have a finite higher energy 
than the lowest energy state within M(Q). 

We first consider a subspace M(d) defined by \ip(a)) = 
A s * (O)IV'(O)) where d has non-zero components for two 
qubits acted on by A s > (0), as shown in Fig. [31 We notice 
that A s * (0) commutes with all terms except two edge 
terms C ei and C e2 . Therefore, solving H ga dget inside 
M.{a) is equivalent to solving 

Af {0)H gad getA s * (0) f = Hgadget + V 

inside A4(0), where V = 2J(C ei + C e . 2 ). 
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FIG. 3. A non-ground-state subspace A4(d). A a * (0) anti- 
commutes with two edge terms C ei and C e2 ■ 



Below, we show that the lowest energy states for 
H gadget + V inside .M(O) have finite higher energy than 
those of Hgadget for appropriate choices of parameters 
[/, t and J. For simplicity of discussion, we neglect a 
constant correction resulting from plaquette term H p by 
writing H gadget = H hop = J2 S h s inside A4(0). Then, one 



may write H' 



gadget 



E, 



#{^ 



/i s + iJ* with 



7J* = 



^ ^ + 2J(C ei +C e2 ) 



S={S* ,S!,S 2 } 

where s\ and s* are connected by e%, and e 2 connects s 2 , 
s* (Fig.©. 

Returning to the A-representation, we note that all par- 
ticles except s*,si,S2 are non-interacting and are gov- 
erned under the same Hamiltonian h s as before. Let us 
denote the lowest energy eigenvalue of h s as E . Noting 
that E is upper bounded by — U, it suffices to show that 
H* > —3U > 3Eq for the existence of an energy gap. 

Let H* =H 1+ H 2 where H, = -t Y<s={s*,s u s 2 }( D \ + 



D.) and H 2 = -t/£ s={s ,, Sl 



M 



0)(r 







2J(C ei +C e2 ). Since one cannot minimize H± and H 2 si- 
multaneously, we obtain a lower bound for H* by finding 
minimal energy eigenvalues for H\ and H 2 individually. 
One can verify that Hi > — 6i by directly finding en- 
ergy eigenvalues of H\. Similarly, one can verify that 
H 2 > min(-3t/ + 4 J, -2U - 4 J). Here, let us choose U 
and J such that J = U/8, and H 2 > ~^U. H' d . has a 
provably higher ground state energy than H ga< ig et when 
H* > -6t- ^ > -3C/ > 3S 0l so we simply set U > 12*. 

This proof may be easily generalized to arbitrary M.(d) 
when U > 16i. 

A drawback of this proof is that a small value of t = 



U/16 gives a weak constant gap for h s and thus the gap 
inside .M(O) is ~ 10 _4 [/. Tighter analysis presented in 
appendix [Bl shows that when J = 0.09C/, t = 0.375C/, the 
system has a quite reasonable energy gap of > 0.075J7 
both inside and outside M(0). 

Particle dimension — In this construction, a gadget 
particle is four-dimensional, and a composite particle is 
eight-dimensional after removing the internal degree of 
freedom for B p . This can be improved through defin- 
ing a similar construction on a triangular lattice, leading 
to six-dimensional gadget particles and four-dimensional 
composite particles. In addition, a more elaborate con- 
struction reduces the dimension of the gadget particles to 
three while keeping the dimension of composite particles 
at four (see appendix |D|) . 

Discussion — Our gadget construction can be general- 
ized to the quantum double model, which may be uni- 
versal for topological quantum computation, in a rather 
straightforward way (see appendix |C|) . We expect that 
similar generalizations are possible for other interesting, 
but highly non-local topologically ordered Hamiltonians. 
In addition, our non-perturbative gadget may find use in 
adiabatic quantum computation and Hamiltonian com- 
plexity problems. 

In our construction, we have heavily taken advantage 
of the fact that the terms being simulated commute. 
Whether non-commuting terms can be simulated in this 
way remains an open question. Perhaps insights from re- 
lated problems in theoretical computer science will prove 
fruitful, opening new connections. 

Conclusion — In this paper, we argue the necessity of 
simulating topological quantum codes non-periurbatively, 
and propose a model which is two-body, exactly solv- 
able, and supports completely localized quasi-particle ex- 
citations. While our construction involves large particle 
dimensions and "gadgcty" interaction terms which put 
it beyond modern experimental capabilities, there has 
recently been remarkable theoretical and experimental 
progress in engineering custom interactions between par- 
ticles |18l - l2C| . We hope our construction will provide a 
stepping stone towards physical realizability of topologi- 
cal quantum codes. 
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Appendix A: Designing the shielding term 



Here, we present an explicit procedure to obtain the shielding term H s hi e id which cancels the edge term H e inside 
the ground state subspace M(0). 
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FIG. 4. Construction of the shielding term, (a) The horizontal edge, (b) A violation of an edge term C e . 

Horizontal edge terms: Let us first consider an edge term C e for a horizontal edge e, connecting two stars s and 
s + x (see Fig. HJa)). We represent C e in the "A-representation" (see Eq. (2) in the main article) inside the ground 
state subspace M(0). Since the edge term C e acts only on |A S ) and |A S+ ^} nontrivially, we will compute the diagonal 
elements (A hor |C e |A /lor ) for \X hor ) = \X S ) <g) \X S+ . X ) where X S1 X S+X = 0, . . . , 7. 

As an example, let us compute the diagonal element for A s = 2 and X s + X = (see Fig. 2Kb)). Then, the edge term 
C e is violated since two qubits on the edge e have different values, and (A/i 0r |C e |Ah or ) = — 1. On the other hand, for 
X s = and A s +s = 0, we have (A/i 0T .|C e |A/j 0r ) = 1 since both qubits on the edge e are in |0). 

By repeating similar analyses for every pair of X s and X s +x , we have: 



(A/i 0r |C7 e |A/j 0r ) — 



A s =2,6 and X s+X =0,4: 

A s ^2,6 and A s+S =0,4 

A s =2,6 and A s+S ^0,4 

A s ^2,6 and A s+S ^0,4 



which we can simplify as: 



(Xhor\C e \X hor ) = (l - 25\ s2 - 2S\ S! e) ■ (2(5 As+£ ,o + 2S Xa+ -A - l) • 

Noting that m s = X s (mod 4) for all s inside .M(O), we see that {Xhor\C e \Xhor) = (Xhor\Ti(m s )T r (m s+ $)\Xhor) within 
M(6) where 

Te(m) = 1 - 25 m> 2, T r (m) = 2S mfi - 1. 

Therefore, a shielding term J ■ Ti(m s )T r (m s+ x) exactly cancels the two-body contribution arising from C e since 
(Aw| - J ■ C e + J ■ Ti(m s )T r (m s+i )\X h or) = inside M{6). 

Vertical edge terms: Next, let us consider an edge term C e for a vertical edge e, connecting two stars s and 
s + y. Then, for |A„ er ) = |A S ) ® \X s +y), one can verify that: 

(A„ er |Ce|A„er) = (1 - 26\ B ,1 - 2S Xs ^) ' (l - 2(5 As + fi ,3 - 2S\ s + ij ^) . 

Then, for a shielding term Td(m s )T u (m s+y ) , we have (X ver \ — J ■ C e + J ■ Td(m s )T u (m s+ y)\X ver ) = inside A4(0). 
Therefore, H e + H shldd = inside M(6). 



Appendix B: Improved bound on energy gap 



In this section we improve the energy gap by improving the energy bound on non- ground-state subspaccs Ai{d). 
Within any subspacc, B p = ±1, and when B p = —1 the energy is raised without affecting other terms. Therefore, 
we need only consider non-ground-state subspaces where B p = 1 and neglect a constant correction from H p . As seen 
previously, solving H gadge t inside A4(d) is equivalent to solving H gadge t + V inside .M(O), where: 




E c e 



e^ er GC 



Here, e contains all edges e : C e \ip(d)) = — 1. We note that within -M((J), for all vertical and horizontal edges: 

C ehor = T e (m s )T r (m s+i ) > T e (m s ) +T r (m s+i ) - 1 
C e „ er . = T d (m s )T u (m s+ y) > T d (m s ) + T u (m a +y) - 1. 

Therefore, a lower bound on H gadge t + V , where 



V'=2j( J2 {Ti(m s )+T r (m s+i )-l)+ J] (T d ( 



m s ) +T u (m s+y ) - 1) 

also serves as a lower bound on H gadge t + V. This is useful as H gadge t + V is one-body in the A-rcprcscntation. 
Organizing terms by stars instead of by edges, we find that: 

Hgadget + V = ^ k„ + ^ ti s , , 

s^e s*£c 

where 



4 = fc a .+2j(af.(2i(m..) - 1/2) + <,(T r (m s .) - 1/2)+ 
2j(af.(T d (m s .) - l/2) + o?.(T u (m..) - 1/2)). 



The coefficients a^ , a r s , , a^, , a^ t =0,1 denote whether C e \ip(d)) = ±\ip(d)) for the left, right, down, and up edges 
of s* . The basic idea of our analysis is to find a lower bound for the energy of h' s * for all 2 4 — 1 = 15 cases where 
of. + a r s , + a*. + a%. > 0. 

We can further tighten analysis by rewriting the above equation: 

h'„ = h s *+2j(a{,{T z {m s *) - 1/2 - fy r ) + a r 8t {T r (m s *) - 1/2 + /3 £r )) + 
2j(a%{T d (m s .) - 1/2 - p du ) + a u s ,{T u {m s ») - 1/2 + /3 du )). 

/3^ r and /3d„ are constants to be optimized. This modification simply redistributes constant energy between stars, 
leaving the total Hamiltonian unchanged; ^ s » a g * = J2 S , a r s , and ^ s , af» = J2 S , a"» . We vary parameters to find an 
optimum at J = 0.09U, t = 0.375U, fa r = 0.25, (i du = 0. For these values, ti 8t > E + 0.25U for all a l 8 , + a r s , + af* + 
a", > 0. Recall that E is the ground state energy of h s . 

Since any non-ground-state subspace M{d) must have at least three stars s* for which this holds, the lowest energy 
of any state in M{d) is at least 0.075C/ above the ground state energy. Likewise, at this value of t the gap to a single 
"vortex" excitation in h s is > 0.0375C/. Since an even number of stars s must be excited in this way, any excited state 
inside A / ((0) must have energy at least 0.0751/ higher than the ground state energy. We combine these two bounds 
to prove a quite reasonable energy gap of at least 0.0751/. 



Appendix C: Gadget for quantum double 

Here, we present a generalization of our gadget construction to the quantum double model. Our construction 
and discussion closely parallel that of the toric code, but are somewhat more complicated. We begin by defining a 
modified version of the quantum double model that we will simulate through a two-body Hamiltonian. Consider an 
arbitrary finite group G (which may be non-abelian) , and consider a qudit with an orthogonal basis {\z} : z G G} 
whose dimensionality is \G\. We define the following group operations: 

L B + \z) = \gz) L a _\z) = Izg- 1 ) 

T*\z) = 5 hit \z) T±\z) = 6 h -i t ,\z) 

Note that, while the group may be non-abelian, L 9 + and Vt_ nevertheless commute, as (gz)h^ 1 = g(zh~ 1 ). We 
consider a system of qudits defined on edges of a square lattice with periodic boundary conditions, where two qudits 
live on each edge in our construction (see Fig. HJa)). The system is governed by the Hamiltonian: 

H = H s + H„ + H e . 
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FIG. 5. (a) The modified quantum double model, (b) A star term, (c) A plaquette term, (d) An edge term. 
The star term is: 

s 3 es 

where A 9 S is represented in Fig. [5jb) and Q is a generating set of group G; any element in G is some product of 
elements in Q. The plaquette term is 

p 

where B p projects onto the subspace where the clockwise product of qudit group elements in a plaquette is the identity, 
as represented in Fig. EJc). The edge term is 



H e — —J / J C e , 



where C e projects onto the subspace where the two qudits on e have the same group element, as represented in 
Fig. [SJd). Again, all terms A S1 B pi C e commute and can be minimized simultaneously. It can similarly be verified 
that the ground states of the original and modified quantum double model can be connected through generalized 
controlled-NOT gates between qudits on edges. Grouping four qudits in each plaquette into a single composite 
particle, A 9 becomes four-body, B p becomes one-body, and C e is two-body Below, we will show how A 9 can similarly 
be simulated through two-body interactions. 

Gadget Hamiltonian: We again add a gadget particle to each star. The gadget particle again has a spin degree 
of freedom m a = 0,1,2,3, but now additionally has a group element degree of freedom: g s g Q where Q is a generating 
set of the group G. Therefore, the gadget particle at position s is described by \m 8 ,g 8 ), and has dimension of i\Q\. 
We replace the four-body star terms A 9 S with two-body interaction terms: 



tlgadge 

H n = - 



H p + H e 



flhop T" ^shield 



-J2_^B p , H e — —J 2 J c e 
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FIG. 6. The hopping term (D 9 ) for the quantum double model. 

The hopping term is Hh op = ^ s h s where 

h s = -U\m s = 0)(m s = 0| - 1 ■ Q s -tJ2(D 9 J + D 9 

see 

( D s) f = J2 \m s + l,g)(m 3 ,g\®A 9 s (m s ) (mod 4) 

m 5 =0,l,2,3 

Qs = ^2 \m s = 0,g s =g)(m s = 0,g s = g'\. 

g.g'&Q 

Terms A 9 (m) are products of two L 9 ± operators as depicted in Fig. [5] Since A 9 (m) are one-body operators when 
plaquettes are viewed as composite particles, hopping terms are two-body The hopping terms effectively induce star 
terms since A 9 = rim=o^s( m )- 

The shielding term is: 

HsMeid = j^ (1 + T e (m s )) (1 +T r (m s+ x)) +8g B g a+ * (1 -I>(m 8 )) (1 ~T r (m s+i )) 1 + 
E ( ( X + T ^ m «)) ( X + T ^s+y)) + S gs9s+ - (1 - T d (m s )) (1 - T u (m s+ y)) ) , 



J 

1 



which will decouple effective interactions to make the model exactly solvable. This shielding term may be derived in 
a similar way as was done for the toric code. 

Decomposition into subspaces: We can analogously decompose this gadget Hamiltonian into subspaces to help 
us solve it. We denote computational basis states whose gadget spin values are all |0)s as 

\lp(g,d)) = (^) \m s = 0,9s) gadget <8> \d}qudit- 

S 

We define the subspace M(d) such that it is spanned by all the states which can be reached from \ifj(g, d)} by applying 



10 

(D 9s ) for some g: 

M(d) = / J] ((D 9 /)^ H'(gJ)), for all g, A Y 

and can verify that yVf (d) is an invariant subspace of H ga d ge t ■ Then, we can solve the gadget Hamiltonian inside each 
subspace A4(d) independently. 

Ground state subspace: We solve for the lowest energy state inside .M(O) where the identity element in the 
group G is denoted by 0. We note that B p — 1 inside .M(O) and thus need not be considered. Denoting the total 
number of stars as N, we may view M(0) as the Hilbert space of N particles, using a somewhat more complicated 
"A-representation" : 

(g) \x s ,9s, is) = n ((^ a ) f ) Xs m- W& o))- 

s s 

We note that this representation is redundant by seeing: 

{(D^) A \^§,d)) = A^(g,d)) 
Translating into the A-representation, the above redundancy becomes: 

\\ s + 4,g s J s ) = \\ s ,g Sl g s -f s ). (CI) 

This allows us limit ourselves to A s = 0, 1, 2, 3 in the A-representation, giving each "particle" a Hilbert space of finite 
dimension 4|C?||G|. 

We can confirm that, within M(0), H s hieid + H e =0. Therefore, in the "A-representation", H ga d g et acts as a 
one-body Hamiltonian H gac i get = const + ^ h s , where: 

h„ = -UJ2 \X s =0,9s,f s )(K=0,g s ,f s \ 

f " ' 

~t ^2 E l As + 1 '9sJs)(^s,gs,fs\ +h.c. 

A s =0,l,2,3 g„ U {DS .» +D >. ' 

-t Y, \K = 0,g a ,f a )(X a = 0,g' s ,f s \, 



g„ si, h Qs 

and Eq. (|C1[) is implicit. 

We can write the lowest energy state inside M(0) as \ipcs(®)) = ® s \ a o), where \a ) = J2x g / a o(A)|A, g, /}, and 
returning from the A-representation, we write it as: 

g s f \=0 

One can verify that \ipGs{^)) is the ground state of the gadget Hamiltonian, and the energy gap may be proven in 
the same manner as was done for the toric code. 

Unitary Connection: The ground state \ipGs(0)) is connected to the ground state of the modified quantum 
double model through the following local unitary transformation: U = J\ U s where: 

3 

U S =J2 J2 \ m ^9s)(m s ,gs\ n Af(m) f . 

QsZlG m s — m<m s 

In particular, we have U\ip GS (0)) = \aa)®a dget ® \ipQD(0)) qu dit where \d ) = J2 g eg Em=o a o( m )\ m ,g), and 
Vl^QD^Ylqudit = lis (S/eG^s ) |0) wn i cn is a ground state of the modified quantum double. In fact, U maps 
each ground state of the gadget Hamiltonian to a corresponding ground state of the modified quantum double. 
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Particle dimension: The particle dimension can be somewhat reduced in the same way as in the toric code; 
defining a similar construction on a triangular lattice and removing the degree of freedom B p gives us plaqucttes with 
a |G| 2 -dimcnsional Hilbert space and gadget particles with a 6 1<5| -dimensional Hilbcrt space. A modification similar 
to the one presented in appendix IDl can reduce the dimension of the gadget particles to 1 + 2\Q\ while keeping the 
dimension of the composite particles the same. 

Appendix D: Reducing particle dimension 

Here we present an extension of our gadget construction which has three-dimensional gadget particles and four- 
dimensional composite particles. Our construction and discussion closely parallel those of the main paper, but are 
somewhat more complicated. We define a modified version of the toric code on a triangular lattice with periodic 
boundary conditions as represented in Fig. [3 governed by the Hamiltonian: 



s p e 

j£s j£p j£e 



(b) A s (c) . B p 

\ x. x - x .x A 

4 :mi ,f+* 




• • (d) Z^e 

/x x\ o z 

FIG. 7. (a) Construction of the modified toric code on a triangular lattice, (b) A star term, (c) A plaquette term, (d) An edge 
term. 

We group the three qubits in each plaquette into a single composite particle with an eight-dimensional space. While 
B p becomes one-body, and C e is two-body through this grouping, the star term A s is only reduced to six-body, and 
must be simulated using two-body terms involving gadget particles. 

Gadget Hamiltonian: This time we add six gadget particles to each star s, one at each corner i = . . . 5. Each 
gadget particle st has three possible spin values, m Si — 0, 1, 2. We replace the six-body star term A s with two-body 
terms Hhop and H s hi e id which involve the gadget particles: 

Hgadget = H p + H e + Hhop + H shield 

H P = -JJ2 B P> H e = -jJ2Ce- ( D1 ) 

p e 

The hopping term is Hhop = J2s^ s wnere 

h s = -U\m Sl = l)(m Sl =11 — * {D\ +D s )+R-(1- n s f 

5 

D\ =Y^\m Si =2}(m Si =l|0i s (t) + K =0)(m Bi = 2| ® \m Sz+1 = l)(m Si+1 = 0| 

i=0 

5 

n s = 2^ (1 — l m si = 0)(m Si = 0|) 

i=0 

where U, t, and R are some positive constants, and m Si represents the spin value of the zth gadget particle at star 
s. Terms A s (i) are products of two Pauli X operators as depicted in Fig. [SJb). Since A s (i) are one-body operators 
when qubits in a plaquette are viewed as a composite particle, hopping terms are two-body. We note that n s , the 
number of nonzero gadget particles at star s, is a good quantum number, and by choosing a sufficiently positive value 
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of R, we restrict the ground space and low-energy excitations to the subspace where n s = 1. In this subspace, D\ is 
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unitary and the hopping term will effectively induce star terms A s , since A s = (-Dj) as seen in Fig. [E^b) 



gadget particle 




A.{5) 



FIG. 8. Construction of the hopping term D\. (a) Gadget particles live at corners of stars Si. Labeling of corners i is shown 
(blue online). Qubits live at edges, (b) Terms A a (m) that are tensor products of two Pauli X operators. Each term acts on 
two qubits (red online), depending on spin values of gadget particles. 

The shielding term H s hieid consists of two-body interactions between gadget particles: 

Hshield = J*^2 ( T -uo( s ) T u ( s + wo) +T^u 1 (s)Tu 1 (s + u 1 ) + Tl^ (s)T,a 2 (s + u 2 )) 

8 

T_ fio (s) = 1-2 (\m S2 = 2)(m S2 = 2| + \m, S3 = l)(m S3 = 1|) , T fio (s) = 2 (|m, = l)(m So = 1| + |m 85 = 2)(m S6 = 2|) - 1, 
T_ fll («) = 1 -2(|m Sl = 2)(m Sl = 2| + |m S2 = 1)(?ti S2 = 1|) , T fil (s) = 1 - 2 (\m Si = 2)(m Si = 2| + |m S6 = l)(m S5 = 1|) , 
T_ fl2 (.s) = 1 -2(|m S0 = 2)(m S0 = 2| + |m Sl = l)(m„ = 1|) , T fi2 (s) = 1 -2(|m a , = 2)(m S3 = 2| + |m S4 = l)(m S4 = 1|) , 



where s + Uj are unit translations of a "star" s to its nearest neighbors in the directions uq = x, u\ 

—x+^/3y 



-Vsy 



,u 2 



As we will see below, this choice of the shielding term decouples effective interactions between neighboring 
gadget particles, making the model exactly solvable. 

Decomposition into subspaces: Now, we solve the gadget Hamiltonian in Eq. (|D1[) inside the subspace n s = 1. 
It is convenient to further decompose this space into subspaces. Let us denote computational basis states whose 
gadget particles are at "rest" as: \ip{d)) = \o) gadget ® \d)qubu where: \o) gad get = <8> s (\m So = 1) <g)f =1 \m Sl = 0)). 
Again, \d} qu bit represents spin values \dj) for qubits, and we define the subspace Ai(d) to be spanned by all the states 
which can be reached from \ip(d)) by applying D\: 

M(d) = ( ]\{D\) X ° \Hd)), for all A s >. 



We can verify that M. (d) is an invariant subspace of H ga dget , and solve the gadget Hamiltonian inside each subspace 
M(d) independently. 

Ground state subspace: We solve for the ground state inside .M(O). We note that inside -M(0), B p = 1, and 
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thus plaquette terms need not be considered. Again, we use the A-representation: 

|A>=(g)|A s )=n(^) As l^(0)>. (D2) 

s s 

Noting that (D\) 12 = A s , (-DJ) 24 = /, these particles can be considered to have twcnty-four-dimcnsional Hilbcrt 
spaces, X s = 0, . . . , 23. In the A-represcntation, the hopping term Hh op can be written as a one-body Hamiltonian: 
Hhop = J2s h s where: 

h s = - U(\\ s = 0)(X S = 0| + |A S = 12)(A S = 12|) 

23 



- t Y^ (l A * + 1)(A S | + h.c) (mod 24). 



A s =0 



While edge terms C e are not one-body, inside M(6) they are exactly cancelled by the shielding term; H shie i d + H e = 0. 
This is because C e = T_a . (s)Tu . (m s+ n . ), where j = 0,1,2 for horizontal, upper-right diagonal, and upper-left diagonal 
edges respectively, which connect s and s + Uj . Because of this exact cancellation, the gadget Hamiltonian is one-body 
in the "A-represcntation": H ga d g et ~ const + ^ s h s . 

Because of this, all energy eigenstates inside yVf(O) can be written in the tensor product form \a) = (££) s \a s ) where 
|a a ) = J2\ a s(^s)|A s ). The lowest energy state is IV'Gs(O)) = ® s \&o), therefore, returning from the A-representation, 
we can write the ground state as: 
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I^G5(0))=riE a o(A)( J Dl) A |V'(0)). 



s A=0 

As with the construction in the main section, this model has a finite energy gap, four ground states connected to 
those of the original toric code by local unitary transformations, and supports completely localized quasi-particle 
excitations; all of these may be proved in analogous ways. 

By removing the internal degree of freedom B p from the composite particles, we arrive at an exactly solvable model 
involving only four-dimensional (spin-3/2) composite particles and three-dimensional (spin-1) gadget particles. 



